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The main goal of this paper is an appHcation of Bayesian model com- 
parison, based on the posterior probabilities and posterior odds ratios, in 
testing the explanatory power of a set of competing GARCH (ang. Gener- 
aHzed Autoregressive Conditionally Heteroscedastic) specifications, all with 
asymmetric and heavy tailed conditional distributions. In building compet- 
ing volatility models we consider, as an initial specification, conditionally 
Student-t GARCH process with unknown degrees of freedom parameter, 
proposed in [Jj- By introducing skewness into Student-t family and in- 
corporating the resulting class as a conditional distribution we generated 
various GARCH models, which compete in explaining possible asymme- 
try of both conditional and unconditional distribution of financial data. 
In order to make Student-t family skewed we consider various alternative 
mechanisms recently proposed in the literature. In particular, we apply the 
hidden truncation mechanism (see P|, ^), an approach based on the in- 
verse scale factors in the positive and the negative orthant (see ^01); order 
statistics concept (511)j Beta distribution transformation (QS]) and Bern- 
stein density transformation (see [2HI)- Additionally, we consider GARCH 
process with conditional a-Stable distribution, see [HOI, HHI- Based on 
the daily returns of hypothetical financial time series, we discuss the re- 
sults of Bayesian comparison of alternative skewing mechanisms applied in 
the initial Student-t GARCH framework. Additionally, we present formal 
Bayesian inference about conditional asymmetry of the distribution of the 
daily returns in all competing specifications on the basis of the skewness 
measure defined in 0. 
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1. Introduction 



The presence of both, conditional and unconditional skewness (asymme- 
try) of the distributions of the financial time series returns has been rec- 
ognized for decades. But, as suggest pH], only a few attempts to specify 
formally this feature have been made. A proper modeling of skewness in 
the distribution of financial returns is important at least for two reasons. 
Firstly, uncaptured skewness clearly affects inference about all parameters 
of the sampling model. As a consequence the final conclusions, drawn from 
the sampling model which does not allow for asymmetry can be misleading. 
Lanne and Saikkonnen present in the impact of the conditional skew- 
ness assumption on the results of making inference about the volatility and 
expected return. They presented empirical analysis, which showed, that a 
positive and significant relation between return and risk can be uncovered, 
once an appropriate probability distribution is employed to allow for condi- 
tional asymmetry. Motivating the importance of asset pricing model that 
incorporates conditional asymmetry, ^H] emphasize that systematic skew- 
ness is economically important and commands a risk premium. Investigating 
the influence of the assumption of asymmetric distributions in portfolio se- 
lection, [18j concluded, that, if investors prefer right-skewed portfolios, then 
for equal variance one should expect a "skew premium" to reward investors 
willing to invest in left-skewed portfolios. Secondly, in pricing the deriva- 
tives and in risk management, the accurate models, which describe the re- 
turn process are particularly desired. The importance of the assumption of 
conditional skewness in models used for option pricing was presented p?T] . 
[T7| and 0. Additionally, conditional skewness clearly influences the results 
of risk assessment built on the basis of the Value at Risk (VaR) concept. 
Application of time varying volatility models with conditional asymmetric 
distributions in Value at Risk prediction present for a Bayesian approach 
to VaR calculation see p9j . 

Within GARCH (Generalised Autoregressive Conditionally Heteroscedas- 
tic) framework, initially proposed by as a conditionally normal stochastic 
process, fat tailed and possibly asymmetric distributions have been also 
proposed and applied. Osiewalski and Pipieh in [2^1 defined GARCH pro- 
cess with conditional skewed Student-t distribution, which is an asymmet- 
ric generalisation of Student-t family proposed by PHj- In |2^, (2111 and 
[301 GARCH process with conditional a-Stable distribution was considered. 
Some other processes with asymmetric conditional distribution were applied 
in ^ni; [nOj; |n2j; and jH). Despite of the fact, that many researchers found 
the conditionally skewed volatility models better than those, which do not 
allow for asymmetry, there is very hard to find the result of the formal 
comparison of explanatory power of such specifications. Many authors con- 
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elude the superiority of conditional skewed models on the basis either of 
the asymptotically based statistical significance of the skewness excess (see 
e.g. |n2]; or Pni) or of informal likelihood inference (see e.g. (213, IS21)- 
Hence more formal approach to investigating the explanatory power of con- 
ditionally skewed models seems to be necessary. Additionally, the results of 
formal comparison of competing unnested specifications of conditional skew- 
ness could be very valuable in selection the best skewing mechanism. 
On the other hand, in recent years in statistics it can be noticed a peculiar 
interest in the theory and applications of distributions that can account for 
skewness. This resurgent field of research yields new families of possibly 
asymmetric sampling models, as well as more general methods of measur- 
ing skewness phenomenon. The most common approach to the creation 
of the family of skewed distributions is to introduce skewness into an origi- 
nally symmetric family of distributions. This approach underlies the general 
classes of skewed probability distributions generated for example by hid- 
den truncation mechanism (see (2], 0), inverse scale factors applied to the 
positive and the negative orthant (see ^Hl), order statistics concept (|14|'l. 
Beta distribution transformation (^3), Bernstein density transformation 
(see [2H]) and the constructive method recently proposed by 
The main goal of this paper is to define a set of competing GARCH spec- 
ifications, all with asymmetric conditional distributions, which also allow 
for heavy tails. As an initial specification we consider GARCH model with 
conditional Student-t distribution with unknown degrees of freedom param- 
eter, proposed by By introducing skewness, according to the methods 
mentioned above, and by incorporating the resulting family as a conditional 
distribution, we generate GARCH models which compete in explaining pos- 
sible asymmetry of the conditional and unconditional distribution of the 
financial data. We also consider GARCH process with conditional a-Stable 
distribution, which, from the definition, also allows for skewness, see |25| . 
By application of Bayesian approach to model comparison, based on the pos- 
terior probabilities and posteriori odds ratios, we test formally the explana- 
tory power of competing, conditionally fat tailed and asymmetric GARCH 
processes. Based on the daily returns of hypothetical financial time series, 
we discuss the results of Bayesian comparison of alternative skewing mech- 
anisms and also check the sensitivity of model ranking with respect to the 
changes in prior distribution of model specific parameters. Additionally we 
present formal Bayesian inference about conditional asymmetry in all com- 
peting specifications on the basis of the skewness measure defined in (2]. 
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2. Creating asymmetric distributions 

Let consider parametric family of absolute continuous real random vari- 
ables I = {£f',£f ■ ^ — > K}, parameterized by the vector 9. For each 
value of G ©, by f{-\0) and F{.\6) we denote the density and cumulative 
distribution function (cdf) of £/. Let assume, that for each 9 G Q the den- 
sity f{.\9) is unimodal and symmetric around the mode. Consider another 
parametric family P of absolute continuous random variables, which distri- 
butions are defined over the unit interval, P = {e/;e/ : ^ (0,1)}, with 
density p{-\r]p) parameterized by vector r]p G H. The unified representation 
of univariate skewed distributions that we study in this paper is based on 
the inverse probability transformation. In our approach the class / is the 
initial family of symmetric distributions, while the class P defines formally 
skewing mechanism. The family of absolute continuous random variables 
IP = {es,£s ■ ^ — *■ I^}, with general form of density s{.\9,rip) is said to be 
the skewed version of the symmetric family /, if the density s is given by 
the form: 

s{x\9,7]p) = f{x\9)-p{F{x\9)\r]p), for xeR (1) 

A number of simple but very powerful results can be obtained from decom- 
position see The most important and rather intuitive fact is that 
the distributions s and / are identical if and only if p{.\r]p) is the density of 
the uniform distribution over the unit interval; i.e. if p{y\rip) = 1, for each 
y £ (0,1). Hence if we want to create the family of distributions IP such 
that I C IP, we must assure, that the uniform distribution over (0, 1) can 
be obtained in family P for some specific value rj* £ H . 
Within the general form several classes of distributions P have been con- 
sidered and incorporated into some specific families of symmetric random 
variables in order to obtain skewness. The first approach of making dis- 
tribution F{.\9) skewed applied hidden truncation ideas. The skew-Normal 
distribution in |^ constitutes the first explicit formulation of such a mecha- 
nism. In general this approach assumes, that: 

s{x\9,-i2) = 2- f{x\9)F{-i2-x\9), for xGM (2) 

where 72 G M is the only one parameter which governs the skewing mecha- 
nism; Tjp = (72). In this case, it can be shown, thatp(?/|72) = 2F(72F~^(y)|0), 
for y G (0, 1). In Q positive and negative values of 72 define right and left 
skewed distributions. Since, for each y G (0,1), it is true that p(y|0) = 
2F(0F~-'^(y)|72) = 1, the case 72 = retrieves symmetry. As an alternative 
it was proposed in [14J to apply the family of Beta distributions in order to 
define p(.|r/p). In particular, s(x|^,73) can be defined as follows: 

s(3;|^,73) = /(x|e)i?e(F(x|0)|73,73"'), for xGM (3) 
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where-Be(y|a, h) is the value of the density function of the Beta distribution 
with parameters a > and 6 > 0, calculated at y E (0, 1). Since i3e(.|l, 1) 
defines the density of the uniform distribution, we obtain, that for 73 = 
7;^^ = 1 the density s is symmetric. In © there is still only one parameter 
73 > 0, which defines the type of asymmetry. If 73 > 1, then s is right 
asymmetric, while 73 < 1 constitutes left asymmetric density. 
The family IP of skewed distributions proposed in (jSJ can be generalized, by 
imposing Beta distribution transformation with two free parameters a > 
and 6 > 0. This leads to the following form for s 

s{x\e,7]p) = f{x\e)Be{F{x\9)\a,b), for xeR (4) 

In this case the vector rjp = (a, b) contains two parameters, which govern 
skewness. As a consequence such a mechanism enables to vary tail weight. If 
a = 6 = 1 we go back to symmetry, while a < b or a > b defines left or right 
skewness. It can be shown that the skewing mechanism in case when 
/ is the family of Student-t distributions, yields skewed Student-t family of 
distributions proposed in advance in |15| . 

Another method for introducing skewness into an unimodal distribution is 
based on the inverse scale factors on the left and on the right side of the 
mode of the density f{-\0). Investigating this concept Fernandez and Steel 
proposed in (H] skewed Student-t family of distributions with the density 
0, 1, 71) defined as follows: 

2 

fsksix\i^,0, 1,71) = —{Mx-fi\u, l,0)/(_oo,o)+/t(a;7rV>liO)/(o,+oo)} 

71 +7i 

where ft{z\h', 1,0) is the value of the density function of the Student-t dis- 
tribution with f degrees of freedom, zero mode and unit inverse precision, 
calculated at z S M. The approach studied in |JU^ can be applied to any 
family I of symmetric distributions by defining in the following skewing 
mechanism for each y G (0, 1): 

,1 , 2 {/hlF-'(9))/(0,(,.5)+/(7r'F-'(!/))V&l)) „, 

"("1^') = ^TT^ nf^) ■ 

where 71 > 0. The resulting density s(.|0,7i) is symmetric if 71 = 1, while 
71 > 1 or 71 < 1 make distribution right or left skewed. 
As pointed in p2] the general form of density s in Q seems to be the 
good starting point in completely nonparametric treatment of the skewing 
mechanism p. As : — > (0, 1) can be in general any random variable 
with probability distribution defined over the unit interval, the possibility 
to model it in an unrestricted fashion is tempting. The next approach of 
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constructing p is a compromise between totally flexible skewing mechanism 
and one obtained in parametric fashion. It uses Bernstein densities (see 
e.g. |2Hl)i which are convex discrete mixtures of appropriate densities of 
Beta distributions. The following form on p constitutes another skewing 
mechanism: 

m 

p{y\wi, Wm) = ^ WjBe{y\j, m-j + 1) , y 6 (0, 1) 
i=i 

where m > 0, Wj > 0, wi + . . . + Wm = 1- The resulting s{.\6,rip) takes the 
form: 

m 

s{x\9,r^p)= f{x\e)-Y,WiBe{F{x\e)\j,m-j + l) /or x E M (6) 
i=i 

where r/p = {wi, . . . , Wm-i), Wj G (0, 1) for j = 1, . . . , m — 1, and in © 
Wm = 1 — wi — . . . — Wm-i- For any m > 0, if Wj = m~^, for each 
j = — 1, then Bernstein density reduces to the uniform distribution 

case. Hence equal weights Wj lead to the symmetry in 

In the next section we present basic model framework, which is a starting 
point in generating conditionally heteroscedastic models for daily returns. 
In order to create the set of competing specifications, we make use of all 
presented skewing mechanisms. We also consider GARCH process with con- 
ditional a-Stable distribution. 



3. Basic model framework and competing skewed conditional 

distributions 

Let denote by Xj the value of a currency at time j. Following f5J,|^ , 
[2ZI let consider an AR(2) process for \n.Xj with asymmetric GARCH(1, 1) 
error. In terms of logarithmic growth rates yj = 1001n(xj/xj_i) our basic 
model framework is defined by the following equation: 

Vj - 5 = p{yj-i- 5) + h^'a.Xj^i + Ej j = l,2,... (7) 

The AR(2) formulation adopted from enables to make inference on the 
presence of a unit root in Inxj. If 5i = 0, then lO reduces to the AR(1) 
process for y^-, i.e. an 1(1) process for Inxj. In an initial specification Mq 
we assume, that the error term Sj = Zj{hj)^'^, where zj are independent, 
Student-t random variables, with u > degrees of freedom parameter, mode 
Ci £ (— oOj+oo), and unit inverse precision; i.e. Zj ~ iiSt{i'XiA)- The 
density of the distribution of the random variable Zj is given as follows: 

-(i/+l)/2 

(8) 
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Defining hj we follow GJR-GARCH(1, 1) specification proposed in jT2]: 

hj = aQ+aie]^J{ej-i < 0)+aie]_iI{ej.i > 0)+6i/ij_i j = 1,2,... (9) 

which allows to model asymmetric reaction of conditional dispersion mea- 
sure hj to positive and negative sign of shock £j-i. 

As a consequence, in model Mq, the conditional distribution of ej (with 
respect to the whole past of the process, tpj-i = {. . . ,ej-2,£j-i)) is a 
Student-i distribution with > degrees of freedom parameter, mode 
Ci G (— oo,+oo), and inverse precision hj; i.e. ej\tpj-i, Mq ~ iiSt{i',C,i,hj). 
In specification Mq the conditional distribution of yj is the Student-t distri- 
bution with > degrees of freedom parameter, mode fij = 6 + p{yj-i — 
(5) -|- (5i lnxj_i -1- Ci^f^ and inverse precision hj (given by the equation (jHl): 

p{yj\'4'j-i,Mo,e,iy) = ft{yj\fij,hj,iy), j = 1,2,... 

where 9 = {6, p,6i,ao,ai,af ,bi,ho) is the vector of all parameters defined 
in sampling model Mq except the degrees of freedom parameter u. 
Now we want to construct a set of competing GARCH specifications {Mi, i = 
1, . . . ,k} by introducing skewness into conditional distribution of yj in Mq. 
The resulting asymmetric distributions are obtained by skewing the distri- 
bution of the random variable Zj, (jHJ, according to methods presented in the 
previous section. The resulting skewed density of Zj is of the general form 
given by Q: 

p{z\Mi)= ft{z\Q,l,u)p[Ft{z-Ci)\'nuMi], /orZjGM, i = l,2,...,k 

where p{.\rii. Mi) defines the skewing mechanism parameterized by the vector 
rji, and Ft{.) is the cumulative distribution function of the Student-t random 
variable with v > Q degrees of freedom parameter, zero mode and unit inverse 
precision. The resulting conditional distribution of Ej in model Mi takes the 
form: 

p(£,#,_i,M,) = hih-'^-^e, - Ci)|0, l,u)hf-''p[Ft{hf-\e, - Ci))|r/„M,] 

where /t(.|0, l,z^) is defined by the formula ((HI). This leads to the general 
form of the conditional distribution of daily return yj in model Mf. 

p{y,\i;,.i,9,u,rj„Mi) = ft{z*\u,0,l)hj'^-''p[Ft{z*\vr, M,], (10) 

where z* = h~^'^{ej — p-j)- As the first specification, namely Mi, we consider 
GARCH model with skewed Student-i distribution obtained by the method 
proposed in |10]. The skewing mechanism Mi] is given by the formula 

(EJ, where r/i = 71 > , and 71 = 1 defines symmetry (i.e. Mi reduces 
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to the model Mq under restriction 71 = 1). The model M2 is the result of 
skewing conditional distribution p{yj\ipj-i, Mq,6) according to the hidden 
truncation method. In this case p[.\r]2,M2] is defined by ^,t]2 = 72 G I^, 
while 72 = defines symmetric Student-i conditional distribution for yj. 
In model M3 we apply |14| Beta skewing mechanism with one asymmetry 
parameter. The skewing distribution p[.|r/3,M3] is defined by where 
% = 73 > 0, and 73 = 1 reduces our model to the case of Mq. Specification 
M4 is based on the Skewed Student-t distribution proposed by ^H]- In this 
case p[.|774,M4] is defined by the formula 774 = (a, 6), for a > and 
6 > and a = b = 1 reduces M4 to Mq. In model M5 we apply Bernstein 
density based skewing mechanism with m = 2 free parameters. It means 
that the skewing mechanism p[.\r]5,M^] is defined by the formula © and 
V5 = (^1) ^2)- The case wi = W2 = 1/3 defines symmetry of the conditional 
distribution of yj , given M5 . 

As an alternative for all methods of making family of Student-i random 
variables skewed, it is possible to consider in a GARCH framework a class 
of distributions, which directly, from the definition, enables for fat tails and 
skewness. The next GARCH specification is based on the assumption of 
conditional a-stability. In GARCH model Mg, as a specification which is 
not a direct generalization of model Mq, we considered in conditional 
a-Stable distribution. In particular we put ej = Zj{hj)^'^, where zj are 
independent a-Stable random variables with a G (0,2], location parameter 
Ci G (— cxD,+oo), unit scale and skewness parameter (3 G [—1,1]; i.e. Zj ~ 
iiSta{Ci, 1) /3, a). For a report of Bayesian inference in model Mq see |29] . 
We denote by y*^*^ = (yi, . . . ,yt) the vector of observed up to day t (used 

in estimation in day t) daily growth rates and by y^^ = {yt+i, ■ ■ ■ , Vt+n) the 
vector of forecasted observables at time t. The following density represents 
the i-th sampling model {i = 1, 2, 3, 4, 5, 6) at time t: 

t+n 

p{y''^\yf\0,^^ur]i,Mi) = Y[p{yj\'>Pj-i,0,uJi,rji,Mi) i = 1, . . . ,6, 

i=i 

where uji is the vector of additional parameters of the sampling model, which 
are not included in 6 and rji; for each i = 1,2, 3, 4, 5 = v, while ujq = a. 
The sampling model Mj is based on the product of the appropriate densities 
p{yj\ij:j-i,9i,uji. Mi), which are generally specified in the formula ifTT)]) for 
i = 1,2,3,4,5, while in case i = 6 p{yj\'ipj-i,9,iVQ,r]Q, Mq) is defined by the 
appropriate density of a-Stable distribution (see |29|). 
Constructed at time t Bayesian model Mi, i.e. the joint distribution of the 
observables {y^^\yj^) and the vector of parameters {6,u>i,r]i) takes the form: 

p{y^'\yf,e, uJi,m\Mi) = p{y^'\yf \e, uji, r]i,Mi)pi9, u;i,7ji\Mi) (11) 
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and requires formulation of the prior distribution p{9^uJi,rji\Mi), for each 
specification Mj, for i = 1,2,3,4,5,6. In general we assumed the following 
prior independence: 

p{e,oj^,7^i\Mi)=p{e\Mi)p{uj^\Mi)p{r]i\Mi) i = l,2,...,6 (12) 

The prior information about the common parameters 9 was initially for- 
mulated by [27.\. For i = 1,2,3,4,5 the prior density p{oJi\Mi) = p{v\Mi) 
defines exponential distribution with mean 10 for the degrees of freedom 
parameter v . In case of conditionally a-Stable GARCH model {i = 6) 
the density p{ijJq\Mq) = p{a\MQ) defines the uniform prior distribution over 
the interval (0,2] for the index of stability a. For i = 1, ryi = 71 > 0, and 
p(ryi|Mi) is the density of the standardized lognormal distribution truncated 
to the interval 71 S (0.5; 2). For z = 2, 7/2 = 72 £ IR, and p{rj2\M2) is the 
density of the normal distribution with zero mean and variance equal to 3. 
For i = 3, 1^3 = 73 > 0, and p{r}z\M'i) is the density of the standardized 
lognormal distribution. In case of i = 4, 1^4 = (a, 6), and p(??4|M4) is the 
product of the densities of the standardized lognormal distribution. For 
i = 5, 7?5 = {wi,W2) and p{ri^\M^) is the product of the normal densities, 
both with mean 0.33 and variance 36, truncated by the following set of re- 
strictions: wi > 0,it)2 > 0, +W2<'\-. For i = 6 r]Q = P , and p{r]Q\MQ) is 
the density of the uniform distribution over the interval [—1, 1]. 

4. Empirical results 

In this part we present an empirical example of Bayesian comparison of 
all competing specifications. As a basic dataset we considered T = 1398 
observations of daily growth rates, yj, of the WIBOR one month Zloty in- 
terest rate from 20.03.97 till 05.09.02. The variability of daily returns i/j as 
well as some descriptive statistics are presented on Figure 1. It is clear, that 
dynamics of daily returns of the WIBORlm instrument is very anomalous. 
Huge outliers, caused by changes in the monetary policy, together with the 
regions of almost no variability, depicts very volatile behavior of rates of daily 
changes of the Polish zloty middle term interest rate. In spite of the fact, 
that in five years from March 1997 to September 2002, the Polish money 
market was changing, our first attempt to compare all models was based 
on the whole dataset. As seen in Figure 1, negative value of the skewness 
statistics clearly shows substantial asymmetry of the empirical distribution. 
It also may indicate skewness of the conditional distribution of yj . 
In Table 1 we present the results of Bayesian comparison of explanatory 
power of all competing specifications. In rows we put the decimal logarithm 
of marginal data densities p(y^*^|Mj) (i = 0, 1, . . . , 6), posterior probabili- 
ties of all models including Mq, posterior probabilities of all conditionally 
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asymmetric GARCH specifications (Mj, i = 1, . . . ,6) and Bayes factors of 
Mq (representing conditional symmetry) against Mj, i = 1,...,6 (repre- 
senting alternative i.e. conditional asymmetry). Both sets of P(Mj|y(*)), for 

1 = 0, 1, . . . , 6 and i = 1, . . . , 6 were obtained by imposing equal prior model 
probabilities. 

It is clear, that the modeled dataset of daily returns of WIBORlm inter- 
est rates do not support decisively superiority of any of competing skewing 
mechanism. The mass of posterior probabilities is rather dispersed among 
models. However, the greatest value if P(Mj|y(*)) receives conditionally 
skewed Student-t GARCH model built on the basis of the hidden trunca- 
tion idea. In this case the value of posterior probability is greater than 
44%. The considered dataset also support conditionally skewed Student-t 
GARCH model with Beta distribution transformation (M3) and condition- 
ally a-Stable GARCH specification {Mq). Those three models cumulate 
about 85% of the posterior probability mass, making all remained speci- 
fications (including conditionally symmetric Mq) rather improbable in the 
view of the data. Very small value of posterior probability is received by 
model M4, which, just like M3, is built on the basis of the Beta distribu- 
tion transformation, but with two free parameters governing the type of 
skewness. The observed time series support parsimony of Beta distribution 
transformation with one free skewness parameter (in M3) and rejects gen- 
eralization proposed by ^Hj- Finally, model M4 receives less than 9% of 
posterior probability mass. Also the Bernstein density transformation (with 

2 free parameters) leads to very doubtful explanatory power of the resulting 
GARCH specification. The model M5 is strongly rejected by the data, as the 
value of posterior probability is more than 10 times smaller than posterior 
probability of symmetric GARCH model (Mq). This may lead to the similar 
conclusion, as it was pointed by |TT], that Bernstein densities do not yield 
flexible skewing mechanism for small values of m, see 

On the basis of posterior odds ratios B^i (for i = 1,...,6) we carried out 
Bayesian testing of conditional asymmetry within presented GARCH frame- 
work, according to the Jeffreys scale, see |^. Except for M5, posterior odds 
Pqi reject conditional symmetry in favor of skewness of the conditional dis- 
tribution of Uj in model Mj. In case of model M2 and M3, the data strongly 
support conditional asymmetry, because P02 and Pqs reaches the third grade 
of Jeffreys scale. The data substantially (grade 2) support M4 and weakly 
(grade 1) support Mi both against symmetric Mq. Additionally, poor ex- 
planatory power of specification M5 is confirmed. The data strongly support 
(grade 3) symmetry against skewing mechanism built on the basis of Bern- 
stein density transformation. 

In Table 2 we present the results of Bayesian inference about tails and skew- 
ness of the conditional distribution of daily returns in all competing speci- 
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fications. The tails oi p{yj\il)j-i^9,u!i,rji^Mi) is modeled by the degrees of 
freedom parameter > in Mj, for i = 0, 1, ...,5, while for z = 6 they are 
captured by the index of stability a G (0,2]. Apart from making inference 
about model specific skewness parameters in all models, we also put poste- 
rior means and standard deviations of skewness measure 7m, proposed by 
[2]. We also present values of posterior probability of left asymmetry of the 
density p{yj\i{jj-i,6,uji,r]i, Mi) (i.e. P(7M <0|y(*),Mi)). 
In case of conditional symmetry (model Mq) the dataset clearly rejects the 
hypothesis of existence of the variance of the distribution p{yj\'ipj-i, 9, u, Mq), 
because the whole density of the posterior distribution of the degrees of 
freedom parameter is located on the left side of the value = 2. Also, 
very tight location of p{u\y^^\Mo) around the value = 1.55 in model 
Mq, assures that the conditional distribution of daily returns possesses first 
moment. Those properties of the posterior distribution p{u\y^^\MQ) re- 
mains practically unchanged in case of the majority of conditionally skewed 
specifications. Only in model M4, Beta distribution transformation with 2 
free parameters substantially changes both, location and scale of the poste- 
rior density of the degrees of freedom parameter. In spite of the fact, that 
p{u\y^^\M4) is located on the right side of the value z/ = 2, the posterior 
standard deviation (equal to 0.55) leaves great uncertainty about existence 
of the second moment of the conditional distribution p{yj\ipj-i6,iy,r]4, M4). 
The similar conclusions can be drawn given model Mg, i.e. given the con- 
ditionally a-Stable GARCH specification. Since the posterior mean of the 
index of stability a locates p{a\y^^\MQ) around the value a = 1.21 (with 
posterior standard deviation 0.04), the dataset decisively rejects conditional 
normality in model Mq (corresponding to a = 2). Prom the definition of 
the family of the a-Stable distributions the resulting conditional distribution 
p{yj\ipj-i, 9, a, r/6, Mq) does not have variance (just like in Mj, i = 0,1, 4). 
Also, posterior distribution of a is located on the right side of the value 
a = 1. It clearly assures the existence of conditional mean of the distri- 
bution of modeled daily returns (again just like in Mi, i = 0,1, ...,5). The 
posterior means and standard deviations of both, asymmetry parameters rji 
and skewness measure jm indicate, that in all specifications Mi, i = 1, . . . ,6 
there is quite strong evidence in favor of left skewness of the conditional 
distribution of modeled daily returns. The posterior distributions of 7m are 
located on the left side of the value 7m = 0, confirming left asymmetry of 
p{yj\'4}j-i,9,uJi,rii,Mi). However, relatively great values of posterior stan- 
dard deviations of 7m reduces potential strength of conditional skewness ef- 
fect. As measured by posterior mean oip{'^M\y^^\ Mi), the greatest intensifi- 
cation of skewness olp{yj\ijjj-i,6, oJi,r]i, Mi) is obtained in model Mi. In this 
case of GARCH model the posterior expectation of asymmetry measure is 
equal to M = —0.063, with posterior standard deviation equal about 0.033. 
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All remained conditionally skewed GARCH specifications generated poste- 
rior distributions of 7Af , localized much closer to the value 7m = and also 
much more dispersed. As a consequence, model Mi yields the greatest value 
of posterior probability of left asymmetry of p{yj\ipj-i,6, iOi,r]i, Mi). In case 
of models Mi, M2, M3, Mg the posterior probabilities P{jm < Oly*-*-*, Mj) are 
greater than 91%. The conditionally skewed GARCH specification based on 
the Berstein density transformation (M5) generates relatively low value of 
posterior probability of left asymmetry, making symmetry, as well as skew- 
ness to the right not strongly rejected by the data. Given Beta distribution 
transformation with two free parameters, posterior probability of 74/ < is 
much lower than in case of model M3. Again we may conclude, that gen- 
eralization, based on two free parameters in Beta distribution, substantially 
changes inference about the properties of p{yj\ipj-i9,uJi,r]i, Mi). Finally, 
on the basis of the Bayesian model pooling technique, we obtained poste- 
rior probability of left asymmetry calculated considering the whole class of 
specifications Mj, i = 1,...,6. The modeled dataset clearly supports left 
asymmetry, as P{'^m < Oly*-*-*) = 0.9263, but it also leaves some uncertainty 
about true intensification of this phenomenon. Posterior probability of sym- 
metry and right skewness of p{yj\ipj-i,9,u;i,r]i, Mi) (equal to 0.0737) does 
not totally reject those cases. 



5. Concluding remarks 

The main goal of this paper was an application of Bayesian model com- 
parison in testing the explanatory power of the set of competing GARCH 
(ang. Generalized Autoregressive Conditionally Heteroscedastic) specifica- 
tions, all with asymmetric and heavy tailed conditional distributions. As an 
initial specification we considered GARCH process with conditional Student- 
t distribution with unknown degrees of freedom parameter, proposed by [Zj. 
By introducing skewness into Student-t family and by incorporating the 
resulting class as a conditional distribution we generated various GARCH 
models, which compete in explaining possible asymmetry of both condi- 
tional and unconditional distribution of financial returns. In order to make 
Student-i family skewed we considered various alternative methods recently 
proposed in the literature. In particular, we applied the hidden truncation 
mechanism (see fS^, fT]), an approach based on the inverse scale factors in the 
positive and the negative orthant (see UJij), order statistics concept jT3|, two 
different settings of the Beta distribution transformation ^H] and Bernstein 
density transformation (see f^). Additionally, we presented the results of 
Bayesian inference within GARCH process with conditional a-Stable distri- 
bution, (see 1313, EHI)- 

Analysis of posterior probabilities of competing specifications did not lead 
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to decisive conclusion about superiority of any of the considered specifi- 
cations. The greatest value of P(Mj|y(*)) received conditionally skewed 
Student-i conditional distribution built on the basis of the hidden truncation 
mechanism (see fS^). The data also supported Beta distribution transforma- 
tion with single free parameter and conditionally a-Stable GARCH process. 
Those three models cumulated more than 85% of the posterior probability 
mass. 

The results of Bayesian estimation showed, that in each competing speci- 
fication the modeled data set confirmed left asymmetry of the conditional 
distribution p{yj\7pj-i,6,uJi,r]i, Mi). In all models Mi (i = 1, . . . ,6) the pos- 
terior distribution of skewness measure 7Af was situated on the left side of 
the value 7m = (representing symmetry). However, substantial dispersion 
of p(M|?/(*\ Mj), as measured by the posterior standard deviation of 7m, did 
not preclude symmetry or right skewness of p(yj|V'j-i, 0, Wj, r]i,Mi). As a re- 
sult, the posterior probability of left asymmetry (equal to 0.9263), obtained 
by application of Bayesian model pooling approach, left some uncertainty 
about the true strength of conditional skewness phenomenon. 

References 

[1] Arnold B.C., Beaver R.J. (2002) Skewed Multivariate Models Related to the 
Hidden Truncation and/or Selective Reporting, Test 11, 7-54 

[2] Arnold B.C., Groenveld R.A. (1995) Measuring Skewness with Respect to the 
Mode, The American Statistician 49, 34-38. 

[3] Azzalini A. (1985) A Class of Distributions which Includes the Normal Ones, 
Scandinavian Journal of Statistics 12, 171-178. 

[4] Bauwens L., Lubrano M. (1997) Bayesian option pricing using asymmetric 
GARCH, CORE Discussion Paper No. 9759 Universite Catholique de Louvain, 
Louvain-la-Neuve. 

[5] Bauwens L., Lubrano M., Richard J-F. (1999) Bayesian Inference in Dynamic 
Econometric Models, Oxford University Press, Oxford. 

[6] Bollerslev T. (1986) Generalised Autoregressive Conditional Heteroscedastic- 
ity. Journal of Econometrics 31, 307-327. 

[7] Bollerslev T. (1987) A Conditionally Heteroskedastic Time Series Model for 
Speculative Prices and Rates of Return, The Review of Economics and Statis- 
tics 69, 542-47. 

[8] Christoffersen P., Heston S., Jacobs K. (2003) Option Valuation with Condi- 
tional Skewnes, Working Paper, McGill University. 

[9] Doman M., Doman R. (2004) Ekonometryczne modelowanie dynamiki pol- 
skiego rynku finansowego, Wydawnictwo AE w Poznaniu, Poznari. 
[10] Fernandez C, Steel M.F.J. (1998) On Bayesian Modelling of Fat Tails and 
Skewness, Journal of the American Statistical Association 93, 359-371. 



14 Bayesian.tex printed on 2nd February 2008 



[11] Ferreira J.T.A.S, Steel M.F.J. (2004) A Constructive Representation of Uni- 
variate Skewed Distributions, Department of Statistics University of Warwick 
technical report. 

[12] Glosten L.R., Jagannathan R., Runkle D.E. (1993) On the relation between 
the expected value and the volatility of the nominal excess return on stocks, 
Journal of Finance 48, 1779-1801. 

[13] Jeffreys H. (1961) Theory of Probability, Oxford University Press, Oxford. 

[14] Jones M.C. (2004) Families of Distributions Arising from Distributions of Or- 
der Statistics, Test 13, 1-43. 

[15] Jones M.C, Faddy M.J. (2003) A Skew Extension of the t-Distribution, with 
Applications, Journal of Royal Statistical Association B 65, 159-174. 

[16] Harvey C.R., Siddique A. (2000) Conditional Skewnes in Asset Pricing Models, 
Journal of Finance 55, 1263-1295. 

[17] Kalimipalli M., Sivakumar R. (2003) Goes Skewness Matter? Evidence From 
the Index Options Market, Working Paper, Wilfrid Laurier University. 

[18] Kim T.H., White A. (2004) On a More Robust Estimation of Skewnes and 
Kutrosis: Simulation and AppHcation to the SP500 Index, Finance Research 
Letters 1, 56-70. 

[19] Lanne M., Saikkonen P. (2004) A Skewed G ARCH- in-Mean Model: An Ap- 
plication to U.S. Stock Returns, Unpublished manuscript. 

[20] Liu S., Brorsen B.W. (1995) Maximum Likelihood Estimation of a GARCH 
Stable model. Journal of Applied Econometrics 10, 273-285. 

[21] McCulloch, J.H. (1985) Interest-Risk Sensitive Deposit Insurance Premia: Sta- 
ble ARCH Estimates, Journal of Banking and Finance 9, 137-156. 

[22] Mittnik S., Doganoglu T., Chenyao D. (1999) Computing the Probability Den- 
sity Function of the Stable Paretian Distribution, Mathematical and Computer 
Modelling 29, 235-240. 

[23] Nadarajah S., Kotz S. (2003) Skewed Distributions Generated by the Normal 
kernel. Statistics Probability Letters 65: 269-277. 

[24] Nelson D. (1991) Conditional Heteroskedasticity in Asset Returns: a New 
Approach, Econometrica 59, 347-370. 

[25] Nolan J. (1999) An Algorithm for Evaluating Stable Densities in Zolotarev's 
M Parameterisation, Mathematical and Computer Modelling 29, 229-233. 

[26] Osiewalski J., Pipieh M. (1999) Baycsian Forecasting of Foreign Exchange 
Rates Using GARCH Models with Skewed-t Conditional Distributions, 25-th 
International Conference MACROMODELS'98, ed. Wladyslaw Welfe, Lodz, 
354-369 

[27] Osiewalski J., Pipieh M. (2003) Univariate GARCH processes with asymme- 
tries and GARCH-In-Mean effects: Bayesian analysis and direct option pric- 
ing, Przegl^d Statystyczny 50, 5-29. 

[28] Petrone S., Wasserman L. (2002) Consistency of Bernstein Polynomial Poste- 
riors, Journal of Royal Statistical Association B 64, 79-100. 



Bayesian.tex printed on 2nd February 2008 



15 



[29] Pipien M. (2005) Value at Risk Estimates and Capital Requirements for Mar- 
ket Risk Obtained From GARCH Predictive Densities, Acta Universitatis 
Lodziensis Folia Oeconomica 190, 197-218. 

[30] Rachev S., Mittnik S., (2002) Stable Paretian Models in Finance, J. Wiley, 
New York. 

[31] Theodossiou P. (1999) Financial Data and the Skewed Generalised-t Distribu- 
tion, Management Science 44, 1650-1661. 

[32] Verhoeven P. McAleer M. (2003) Fat Tails and Asymmetry in Financial 
Volatility Models, Faculty of Economics, University of Tokyo, CIRJE-F-211. 

[33] Zolotarev U.M. (1961) On Analytic Properties of Stable Distribution Laws, 
Selected Translations in Mathematical Statistics and Probability 1, 202-211 




18.03.97 23.12.97 02.10.98 16.07.99 26.04.00 12.02.01 21.11.01 06.09.0! 

mcan=-0,07; std. dev.=l,09; skewness=-3,74; kurtosis=48,50 



Figure 1. Daily returns j/j of the WIBOR 1-month Zloty interest rate from 
20.03.1997 to 05.09.2002, T = 1398 observations 
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Table 1. Logarithms of marginal data densities, posterior probabilities of all 
competing models (including Mq) and Bayes factors of Mq against Mj, for i = 
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Table 2. Posterior means and standard deviations of tails and asymmetry 
parameters in all models as well as posterior probability of left skewness of 



